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Abstract 

We investigate the total decay rate of the (ground state) B c meson within the framework 
of the relativistic constituent quark model formulated on the light-front (LF) . satisfactory. 
The exclusive semileptonic (SL) and nonleptonic (NL) beauty and charm decays of the B c 
meson are described through vector and axial hadronic form factors, which are calculated 
in terms of the overlap of the parent and daughter meson LF wave functions. The latters 
are derived via the Hamiltonian LF formalism using as input the update version of the 
ISGW model. The inclusive SL and NL partial rates are calculated within a convolution 
approach inspired by the partonic model and involving the same B c wave function which is 
used for evaluation of the exclusive modes. We predict the partial rates for 74 exclusive SL 
and NL channels and 43 inclusive partial rates corresponding to the underlying b — > c and 
c — > s quark decays. Based on our approach we find T b (B c ) = 0.52 ± 0.02 ps' 1 , T C {B C ) = 
0.98 ±0.07 ps _1 , where the theoretical uncertainty is dominated by the uncertainty in the 
choice of the threshold values at which the hadron continuum starts. For the B c lifetime 
we obtain tb c = 0.63 ±0.02 ps in a good agreement with the prediction obtained using the 
nonrelativistic operator product expansion. We also predict decay rates for many specific 
weak transitions of B c . In particular, for the branching fractions of the B~£ — > J/ipfj^fn, 
B+ — » J/V'vr + and — > J/tp X decays we obtain 1.7%, 0.1% and 13.2%, respectively. 
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1 Introduction 



Weak decays of hadrons containing a heavy quark present the most direct way to determine the 
Cabbibo- Kobayashi-Maskawa (CKM) matrix and to test our present knowledge of the QCD 
confinement scale inside hadrons. Among various heavy flavor hadrons the B c meson, the bound 
state of the be system with open charm and beauty, is particularly interesting. The theoretical 
interest in the study of the B c meson is stimulated by the experimental search at CDF and 
LHC Since last year, ALEPH has a very clean signal — > J/ipfi+v^ with a measured B c mass 
5.96ig;i9 GeV/c 2 |L|. A new preliminary OPAL analysis finds two J/^tt + candidates with 
the masses 6.29 ±0.17 and 6.003 ±0.06 GeV/c 2 . Recently the CDF Collaboration reported the 
observation of B c in 1.8 TeV pp collisions using the CDF detector at Fermilab Tevatron ||. 
The CDF results for the mass and lifetime are M Bc = 6.40 ± 0.39(stat) ± 0.l3(syst) GeV/c 2 
and t Bc = 0AQt° il{stat)±0.03{syst) ps. The physics of B c mesons has stimulated much recent 
works on their properties, weak decays and production cross section on high energy colliders. 
A comprehensive analysis of the be spectroscopy and strong and electromagnetic decays of the 
excited states has been given in M. 

Similarly to D and B mesons, the ground be state is stable against strong or electromagnetic 
decay due to its flavor content and disintegrates only via weak interactions. The weak B c meson 
decays occur mainly through the CKM favored b — > c transitions with c being a spectator, and 
c — > s transitions with b being a spectator. Weak decay properties of the ground state be 
including semileptonic (SL) and various exclusive nonleptonic (NL) modes have already been 
discussed in refs. In refs. ||, [0] the B c lifetime has been estimated on the basis of 

a modified spectator model, where the phase space for the free quark decay is modified to 
account for the physically accessible kinematical region |J, or the b and c quark masses are 
reduced by the binding energy to incorporate bound state effects ||. An estimation of the B c 
lifetime t Bc using a modified spectator model and information gained from the calculation of 
dominant exclusive modes is given in ||. A careful analysis of the B c lifetime, performed using 
the nonrelativistic QCD and including the nonrelativistic corrections up to v 2 /c 2 , has been 
carried in ]TD[ , see also [[□]] . The wide range of lifetimes tb c reported in these papers, reflects 
the uncertainty due to the various model assumptions on the modification of the free decay 
rates due to the bound state effects and the limited knowledge of the heavy quark masses. 

Weak decays of charmed and bottom hadrons are particularly simple in the limit of infinite 
heavy quark mass, where the decay rate of a hadron Hq, containing a heavy quark Q is 
completely determined by the decay rate of the heavy quark itself. In this limit, one might 
expect that T(B C ) ~ T(B°) + T(D°). If this result were hold, the B c lifetime would be rather 
short, namely r Bc w 0.3 ps, and B c decays would be dominated by c — > s decay over b — > c decay 
in the ratio of roughly 4:1. In reality, heavy hadrons are the bound states of heavy quarks 
with light constituents. The inclusion of these soft degrees of freedom generates important 
contributions due to the preasymptotic effects, like the binding effects and the Fermi motion 
of a heavy quark inside the hadron. These effects have a significant impact on the lifetime 
and SL branching ratios. The leading nonperturbative effect is described by a distribution 
function F(x), ("shape function"), which arises as a result of the resummation of the heavy 
quark expansion [12|] and has been also incorporated into phenomenological models of inclusive 
decays, first in refs. fl3f , [fLUR , and lately in refs. []T5|| — [1TS|| . The actual calculation of F(x) is a 



difficult nonperturbative problem which in practice introduces considerable uncertainties in the 
evaluation of the hadron lifetime. An important advantage of the be system is the applicability 
of a quark potential model treatment [|J, M, [I9|. In what follows we assume that, instead of 
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QCD with its complicated dynamics of infinite number of degrees of freedom in the light cloud, 
we consider a constituent bound-state problem of a heavy quark interacting with a lighter one 
via a potential. Then, using the formalism of the LF relativistic quantum mechanics [2~0|, 
it is possible to encode all the nonperurbative QCD effects in a LF quark model wave function 
i[}(x,pj_) of a heavy hadron. The internal motion of a heavy Q-quark inside the heavy flavor 
meson is described by the distribution function \i/j(x,pY)\ 2 , which represents the probability 
to find a heavy quark carrying a LF fraction x = Pq/Ph q of the meson momentum and a 
transverse relative momentum squared p\ = p^_. 

A priory, there is no connection between equal-time (ET) wave function w(k 2 ) of a con- 
stituent quark model and LF wave function if)(x, pj_). The former depends on the center-of-mass 
momentum squared k 2 = |k| 2 , while the latters depends on the LF variables x and p\. However, 
there is a simple operational connection between ET and LF wave functions f22|. The idea is to 



find a mapping between the variables of the wave functions that will turn a normalized solution 
of the ET equation of motion into a normalized solution of the different looking LF equation 
of motion. That will allows us to convert the ET wave function, and all the labor behind it, 
into a usable LF wave function. This procedure amounts to a series of reasonable (but naive) 
guesses about what the solution of a relativistic theory involving confining interactions might 
look like. 

We convert from ET to LF momenta by leaving the transverse momenta unchanged, pj_ = 
kj_, and letting 

1. 
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where Mq = \Jk 2 + rrig + ^jk 2 + m 2 p is the free mass operator with m sp being the mass of the 
quark-spectator. Now we obtain the LF wave function from 
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It can be easily verified that / dx J d 2 p±\ip(x,p 2 L )\ 2 = 1 provided / d 3 k\w(k 2 )\ 2 = 1. It is wave 

o 

functions made kinematically relativistic in this fashion, that were used to calculate the form 
factors of heavy-to-heavy and heavy-to-light exclusive transitions in refs. |20| , . 

A relevant feature of our approach is that both exclusive and inclusive decays are consistently 
treated in terms of the same heavy quark wave function ip{x,p\). So far, this approach has 
been applied only to the exclusive and inclusive partial widths of B° [18], where it has been 



found that the overall picture is quantitatively satisfactory. In this paper we extend previous 
calculations to compute the lifetime and various decay branching fractions of B£. 

The paper is organized as follows. In Section 2, after a general overview of the model, we 
set the framework for our theoretical calculation of heavy meson lifetimes in the LF quark 
constituent model. Here we present general analysis focusing on calculation of inclusive SL 
and NL decay rates. In Section 3 we apply the model to calculate various beauty and charm 
decay rates of the B c meson in parallel with those of the beauty decays of the B° meson. In 
the former case we include both weak annihilation and Pauli interference contributions to the 
total width. Section 4 completes the paper with a summary and conclusions. Technical details 
of the calculation are given in Appendix. 
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2 Description of the model 



Without going into too many details, we present first a short outline of the method of calculation 



of heavy meson partial widths, which is discussed at length in [18| . We introduce the necessary 
definitions and describe the main steps in the calculation of exclusive and inclusive decay widths 
of the heavy mesons. We will make several strong assumptions to obtain a model as simple 
as possible. Yet, we will check that it agrees, within reasonable limits, with available data 
on the SL and NL S°-meson decays, as well as with other theoretical predictions. We thus 
feel confident that this model can be used advantageously to obtain rough estimation for other 
hadrons, such as B c . 

2.1 Kinematics 

Consider the SL decay rates first. Instead of considering the exclusive modes individually we 
will sum (in the representative case of B° decay) over all possible charmed final states X^d 
containing the c-quark. This sum includes hadronic states with a large range of invariant mass 
Mx'- Md < Mx < Mb- For the bottom mesons, the 6-quark mass mj, provides a short-distance 
scale that leads to a large energy release into the intermediate hadronic states. Therefore the 
energy which flows into hadronic system X cd - is typically much larger than the energy scale 
Aqcd which characterizes the strong interactions. It will be valid over almost all of the Dalitz 
plot, failing only in the narrow corner region where the observed mass spectrum is dominated by 
the two narrow D and D* peaks. Accordingly, the total SL rate of the B° meson is represented 
in the following form 

T(B° -> Xcatut) = T(B° -> Diut) + T(B° -> D*iu t ) + T(B° -> X'Jv e ), (3) 

where X' 5d represents the charmed hadron continuum, including also the resonance states higher 
than D and D*. The usefulness of such an expansion rests on large energy release in the inclusive 
decay. 

Contrary to refs. Jl7|], fI5|| , where the exclusive decay rates in Eq. (|3|) were calculated using 
the universal Isgur-Wise (IW) function, our calculations of the exclusive rates are not relied on 
the Heavy Quark Effective Theory (HQET). Instead, they use the hadronic form factors that 
depend on dynamics of specific channels. The relevant formulae are collected in the Appendix, 
see Eqs. flOD , (|Q1) . 



Before calculating the inclusive SL decay rates we briefly recall the necessary kinematics. 
The modulus square of the amplitude summed over the final hadronic states is written as 

\M\ 2 = ^\V QV \ 2 L^W a ,, (4) 

where Vq<q is the relevant CKM matrix element, L af3 is the leptonic tensor 

L a/3 = 8[pe a p Ue i3 + PuePPea ~ 9ap(pePu e ) + ie apjS p]p 5 Ue \, (5) 
with e 0123 = 1, and W a p is the hadronic tensor, 
W a/3 = (2tt) 3 

r n _ n M 

E / &{P* - V - EPi) < B°\jm\n >< n\M0)\B° > ft t^w, (6) 
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Here, Pgo is the 4-momentum of B° and q = p e + p u is the 4-momentum of the lepton pair. 
The hadronic tensor can be decomposed into five different Lorentz covariants: 

W a/3 = {-g a p)Wi + v a vpW-z + ie a/3l su J v 5 W 3 + (v a up + u a vp)W4 + u a upW 5 , (7) 

where v is the 4- velocity of decaying Hq, and u = q/Mu Q - The functions Wi depend on 
two invariants, q 2 and go, where the latter variable is related to Mx by: q = (Mjj q + q 2 — 
M x )/2M Hq . It is convenient to scale all momenta by M Hq , so that q 2 = M^ Q -t, M\ = Mjj -s. 
Then the SL width can be cast into the form 



32 r r \q\ 

T SL = -T \V Q/Q \ 2 J dtm,mlml) J dsj^-G(t, s), (8) 
where the prefactor 



2. 



GlM% 

sets the overal scale of the rate, Q(q 2 ,m 2 ,m 2 ) is the usual triangle function scaled by q 
Q(q 2 ,m\,m\) = \jl — 2A + + Xi., with A± = {m\ ±ml)/q 2 , m li2 being the lepton masses, and 

G(t,s) = 3t(l + \ 1 -2\ 2 )W 1 (t,s)+ | (l + AO^ + ^t ) W 2 {t,s) + 

h 2 t ((1+t- s)W 4 (t, s) + tW 5 (t, s)) , (10) 

with Ai = A+ - 2A^, A 2 = A+ - X 2 _ . In Eq. (|) 

2|q| 



= a(t,s) = J(l+t- s) 2 -At, (11) 



and the limits of integrations in the t — s plane are given by 

mP_Y _ r _ (m 1 + m 2 ) 2 _( M^V 

M H Q M Hq V Hq J 



' mm 



These expressions conclude the kinematical analysis. The next task is the calculation of the 
hadronic structure functions Wi. 



2.2 The LF constituent quark model approximation for W^ v . 

The theoretical treatment of inclusive SL decays of heavy-flavor mesons carries a distinct sim- 



ilarity to deep-inelastic lepton-nucleon scattering - this analogy has been used in refs. [17 



[fL8[ within the framework of the constituent LF quark model. The approach is based on the 
assumption of quark-hadron duality which means that the sum over a sufficient number of 
exclusive hadronic decay modes can be described in terms of partonic degrees of freedom. The 
standard strategy is to represent all the states higher than the ground pseudoscalar and vector 
states by the free-quark approximation starting from some effective continuum threshold 



x 



Then the hadronic tensor W a p is given through the optical theorem by the imaginary part of 
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a quark box diagram describing the forward scattering amplitude in the second order in weak 
interactions []: 

Wop = jf ^ / ^T^PQ'iPqMpQ ~ <?) 2 " m 2 Q >\e(E Q ,)mx,pl)\ 2 , (13) 

where Eqi is the energy of Q' and the tensor lJQ®\ describing the Q — > Q'W transitions, is 
defined analogously to the lepton tensor in Eq. (|]): 

L% Q) = 4{ PQ , aPQ(3 + PQ/f3 p Qa - g af 3(pQ>p Q ) + ie a ^ S PQ>p 5 Q ] • (14) 

In Eq. (|14]) the extra factor 1/2 corresponds to the average over the 6-quark spin projections. 

Equation (|I1J) incorporates all long-range QCD effects in the nonperturbative distribution 

i 

function \ifj(x,p 2 ± )\ 2 , whose normalization is given by J dx J d 2 p±\tp(x,p 2 L )\ 2 = 1, while hard 

o 

gluon corrections can be taken into account by perturbative methods and renormalization- 
group techniques. In Eq. (|T3D the ^-function corresponds to the decay of a Q-quark with 
longitudinal momentum xPh q to a Q'-quark and has two roots in x, viz. 



2 2 , 5(x — x + ) + 5(x — x. 

M H Q \X+- X - 



S[(PQ - qY - mj,] = v A/r2 + i\ ^ | ~\ (15) 

where 



x± = -(1 + t- s) ± 



2 

mi 



2 v-.- .^ T l+i- S ) 2 -« + 4^. (16) 

By the quark masses tuq and trqi we hereafter understand the "constituent" quark masses 
taken from a particular constituent quark model. The root x_ is related to the contribution 
of the Z-graph arising from the negative energy components of the Q'-quark propagator and 
is prohibited by the 0(Eq>) in Eq. (|T3|). We now substitute into the quark tensor L^^, Eq. 
flUD , Pq = xPh q , Pq> = Pq — q and use Eqs. fll3"D and (fl5f) to obtain 

W 1 =F(x+), W 2 = 4 * +F{x+ \ , W 3 = W 4 = -2 l F{x+) W 5 = 0, (17) 

|X-f_ Ju— | X — 

where 

F(x)= / |^(x,p 2 ± )| 2 rfV. (18) 



After substituting Wi from Eqs. ( |17l) into (|T0|) the inclusive SL decay rate, Eq. (|8|), is given 

by 



"h q ,iu = -ToJsl\V qq/ \ 2 J dt$(t) J a(t,s) 



[(1 + \ x )a 2 {t, s) %+ + 3t(l - \ 2 _)]F(x+), (19) 
where we have inserted the factor J$l ~ 0.9 representing the effect of the radiative corrections 



24 



1 In what follows we consider the Q- quark as a virtual particle of the mass rriq = x 2 M'jj Q and neglect the 
transverse momenta of the Q-quark. As a result, the expressions for the SL branching ratios are derived in close 
analogy with deep-inelastic scattering. The same result has been derived [ fl6| using the light-cone dominance 
of SL inclusive i?-meson decays, see also |23]. 
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2.3 Nonleptonic decays. 



The calculation of the NL decay rate closely follows the SL one. We expect Hq decays to multi- 
meson states to proceed predominantely via the formation of a quark-antiquark state, followed 
by the creation of the additional qq pairs from the vacuum. The effective weak Lagrangian, 
e.g. for b — * cuq processes with q — d, s is given by 



G F 

L{n) = -^V cb V uq (ciOi + c 2 2 ), 



(20) 



where 0\ and 2 denote current-current operators with the color-nonsinglet and color-singlet 
structure, respectively: 



O x = (cT»(gI», 2 = (ciT^iqjT^u 1 



(21) 



and r M = (1 — 75)7^- The lepton pair is substituted by a quark pair, and the Wilson coeffients 
q are the perturbative QCD corrections that describe the physics between the W boson mass 
and the characteristic hadronic scale of the process. We shall use the values [25] 



Ci(m b ) = 1.132, c 2 (m fe ) = —0.286; c x (m c ) = 1.351, c 2 (m c ) = —0.631, 



(22) 



obtained at next-to-leading order with the evolution of the running coupling constant being 
done at two-loop level using the normalization a s (mz) = 0.118 ±0.003. The inclusive NL rate 
is given by Eq. (|T9f) with the substitution r Jsl — * 3r r/, where 



V 



2 N C + 1 , 2 



2N r . 



+ c 



with c± = c\ ± c 2 . Following ref. [18[] we let N ( 

1 



2N C 
00 in Eq. 



V 



2 1 2 
c, + c_ 



(23) 

J) in which case 77 is reduced to 

(24) 



For the exclusive two-meson NL decays Hq — ► PP, PV, VV, where P and V are the lowest- 
lying pseudoscalar and vector mesons, respectively, we use the BSW approach |26[, pTfl . There 
are two main ingredients in this approach: 



• one assumes factorization: < MiM 2 | J^J^Hq >=< M 2 | J^|0 >< M±\ J^\Hq > to describe 
H Q -> M X M 2 decay, 

• one employs LF meson wave functions to compute < M\\J^Hq >. 

This method which lacks a firm footing for the B meson, can in fact be justified in the case of 
B c . Deviations from factorization arise from higher Fock components of the B c wave function 
and therefore are of higher order in the nonrelativistic expansion. 

By factorizing matrix elements of 4-quark operators contained in ( p0|) , we distinguish two 
classes of NL decays f27] corresponding to two flavor-flow topologies relevant for our discussion: 



the so-called "tree topology" (class I) and the "color-suppressed tree topology" (class II). The 
class I {external decays) contains those decays where only a charged meson can be generated 
directly from a colour-singlet current. As before, we calculate separately the rates for the 
exclusive two-meson decays of the type B° — > D + n~ . B° — > D + n~ , B° — > D + p~ , B° — > D + p~, 
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etc., and the rates for the decays into multimeson states like B° — > Xj c 7r~, B° — > D + X' nd , or 
5° — > Xj +X^ d . A second class of transitions (internal decays) consists of those decays in which 
the neutral meson is generated from the quark current, like J/ip in the decay B° — > KJ/ijj. 
There is the third class of decays where the external and internal amplitudes interfere. These 
decays play important role in the case of the D mesons, but are less important for B c , and are 
absent for B°. The class I and class II two-meson amplitudes are proportional to the QCD 
coefficients 

ai= Cl (/i) + ^, a 2 = cM + °-j^, (25) 

repectively, with /x = 0(m&). The term proportional to 1/N C arises from the Fierz reordering of 
operators Q 2 to produce quark currents to match the quark content of the hadrons in the initial 
and final state after adopting the factorization assumption. This well-known procedure P7[] 
results in matrix elements with the right flavor quantum number but involve both color-singlet 
and color-octet quantum operators. We shall use the "naive" factorization approximation in 
which one discards the color-octet operators 0. In what follows we neglect the 1/N C corrections 
in Eqs. (ESI) letting a x = C\ and a 2 = c 2 . 

In conclusion, we list the expressions used below to calculate the NL partial widths. The 
exclusive two-meson decay rates are calculated using Eqs. ( |A.12|) -( |A.15| ) given in Appendix. 
The exclusive NL rates for the decays of the type B° — > D + X' ad are given by Eqs. ( |A.5| ),( [A^6D . 
The inclusive NL rate of the type B° — > n~(p~)X dc are given by Eqs. ( |A.18| ),( [A.19|) . Finally, 
the inclusive NL rates, like B° — > X' dc + X' ad , have been calculated using Eq. (|T9"1), where the 
lepton pair is substituted by a quark pair and the prefactor r is replaced by 3T cl or 3r C2 
for the external and internal inclusive decays, respectively. 



3 NUMERICAL RESULTS 

We are now ready to add the various contributions presented above and to estimate the B c life- 
time. We first summarize the input parameters that will be used in our numerical calculations. 
The partial decay rates depend on the following set of parameters: 

• the masses of constituent quarks building up final mesons and multi- hadrons and the 
masses of pseudoscalar and vector mesons; the formers are taken from the ISGW2 model 



13: 



m u = rrid = 0.33, m s = 0.55, m c = 1.82, m;, = 5.2, (26) 
while the latters have been taken from the Particle Data Group publication [p9|1 . Note 



that the constituent quark masses mi, and m c satisfy approximately the relation m;, = 
m c + 3.4 GeV, which is consistent with the known formula relating the pole masses m^ po i e 
and m^poie in HQET; 

the meson decay constants fps and fy of the pseudoscalar and vector mesons to the 
V^-boson. For a pseudoscalar meson P = (qiqz), we define < 0|g27^759i|-P(p) >= ifpsPv 
The decay constant of a vector meson V = (gi^) is defined as < 0| ^27^75^1 |V" >= 
Epjnyfy. The constants fps and fy used in our calculations were taken from ref. |3(| . The 



2 To compensate for the neglected octet operators and other non- factorizing contributions one usually treats 
£ = l/N c in Eq. ( p5| ) as a free parameter which can, in principle, be channel dependent. For a detailed discussion 
seerefs. PI 
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B c meson mass is known with a good accuracy from the current quark model calculations 
calculations, we use M Bc = 6.3GeV; 



• the threshold values at which the hadron continuum starts; in our calculations we 
have adopted two different choices: Mjp = Mp + (case a) and = M v + 

(case b), where M P and My are the masses of the pseudoscalar and vector ground state 
mesons, respectively. 

For the radial wave function appearing in Eq. the Gaussian ansatz of the ISGW2 model 
has been adopted in which the main characteristic is confinement. The oscillator parameters 



have been taken from [ 19[ . Their values along with the values of the meson decay constants 
are collected in Table 1. 



3.1 B° decays 

First, we would like to present an overview of different _B c -decay mechanisms and their relative 
importance as obtained within the framework we are advocating. To this end, we first calculate 
the partial B decay modes corresponding to various underlying quark subprocesses. In Table 
2 we show the B° partial widths for two specific choices of the parameter Mjp . The left set of 
numbers collects numerical results obtained for the threshold value for the hadron continuum 
= Mp + m n} while the right set is given for M^P = My + m n , where in the considered case 
M P = Mp, and My = M D * |. Table 2 also shows the numbers N exc and N imc of (both external 
and internal) exclusive and inclusive channels, respectively, for the b —>■ c decays. Our analysis 
incorporates 54 exclusive SL and NL b — > c decays and 29 inclusive b —>■ c decays including 



two baryon-antibaryon channels. The latters were calculated using the Stech approach |31 
We have also included the CKM suppressed b — > u contributions (with IV^/V^I = 0.08) which 
slightly change the overall results. We take the vector and axial form factors for B° — > D(D*)£vi 



and B° — > ir(p)£i>£ transitions from ref. [TL\, where the fit of the form factors by the functions 

/(g2) = 1 - g 2 /Af -V q*/Ar (27) 

where / = F±, V, A , A±, A 2 have been also presented. The distribution functions Fb(x) 

for the B° meson and Fi(x) for the B c meson are shown in Fig. 1. Note that in the later 

i 

case F c (x) = Fi(l — x). In terms of the mean value < x >= J dx x F(x) and the variance 

o 

i 

a 2 = J dx (x— < x >) 2 F(x), one gets < x >= 0.90, \fa i = 0.06 for B° and < x >= 0.72, 
o 

\fo 2 ~ = 0.09 for B c . The difference in the values of a 2 reflects the difference in the values of 
<P±> = (3h q for B° and B c . 
Numerical values of NL branchings depend upon estimation of QCD corrections which give 
contribution of order of 10% — 20%^. The main uncertainty (~ 10%) is related to the choice of 

3 Note that for the former choice there is, in principle, a danger of double counting, because e.g. the exclusive 
rate B° — * D° accounts partially for the two-particle state Dir, which in our approach is included in the inclusive 
rate 

4 The treatment of radiative corrections in our phenomenological application is associated with large uncer- 
tainties. In a typical hadronic process, there are several mass scales involved: the hadron masses, the quark 
masses, the energy release, etc. Thus, there is an uncertainty in the choice of the "characteristic" scale /i of 
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M ( x ] . For each case we then fix the effective value of \V c b\ by the requirement that the measured 
B° meson lifetime rgo ~ 1.56 ps is obtained. Our aim here is not to establish a new value of 
\V c b\, but rather to illustrate how our approach works. Moreover, imposing the \V c b\ constraint 
strongly reduces the dependence of the predicted value of tb c from the uncertainty related to 
the choice of the continuum threshold Q. 

A few tests of the model form factors may be quoted here: The branching ratios for the 
SL decays B° —>■ D*£u and B° — > Dlv are found to be 5.0% and 1.9%, respectively. The 
corresponding experimental values are (4.64±0.26)% and (1.8±0.4)% p9[| . In the SL B° — > p£v, 
the values for the form factors at q 2 = are V = 0.216, A x = 0.170, A 2 = 0.155. The expected 
branching ratio for this decay is 13.2|V^;,| 2 , which when combined with experimental result 
BR(B° -> piv t = (2.5±8q S ) • 10" 4 H yields \V ub \ = (3.1 ± 0.5) • 10~ 3 . The calculated p/vr ratio 



in exclusive SL B decays is 2.2, as compared with the experimental value of 1.4 ± 0.6. 

Our predictions for the branching ratios of the dominant NL two-body decays of the B° 
meson are given in Table 3. The QCD coefficients a± and a 2 have been left as parameters 
in the expressions for the branching ratios. For comparison, we show the corresponding pre- 



dictions obtained using the Neubert-Stech |30[ and NRSX models |33| and the world average 
experimental results, as recently compiled in [34]. 



It is instructive to compare in details some theoretical predictions for the NL two-body 
decays with the data. Due to the uncertainty in the values of the QCD parameters ai and 
a 2 , we first concentrate on the ratios of branching fractions in which these coefficients cancel. 
From the class I transitions listed in Table 3 we obtain 

BR(B°^D + 7i-) r A 

r ^ bS:dV) ~ omioss - om] - (29) 

Here we use the LF ISGW2 model as our nominal choice, and quote results obtained with 
the models of refs. [30 , [33| in brackets. The experimental results for these quantities R± = 



and 



1.11 ±0.23 and R 2 = 1.15 ±0.34 agree with the predictions and do not distinguish between the 
models because the errors are still large. For the class I transitions to final states which differ 
only in their light mesons we obtain 

*s = 4 ' f g ,'£! - « 103 [1.03, 1.06], (30) 



and 



fl BR(B° -> D* 



j; BR(B° -> D* 



7T 



*< = %' BR(B° W °' 91 [ °- 88 ' °' 88] - (31) 



a process. In principle this is not a problem, since the products of the Wilson coefficients with the hadronic 
matrix elements are scale independent. However, we employ simple model estimations of the matrix elements, 
which do not yield an explicit scale dependence that could compensate that of the Wilson coefficients. Instead, 
the quark model calculations are assumed to be valid on a particular scale, typically /i«l-2 GeV. 

5 In addition to using the measured B° lifetime one could be tempted to use the D° lifetime to eliminate the 
similar uncertainties for the c — > s transitions. However, we refrain from this procedure, since our approach 
applied to charm mesons is less reliable than for B c and probably more qualitative than quantitative. 
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The experimental values are R3 = 0.95 ± 0.24 and R4 = 0.99 ± 0.250, again in agreement with 
our predictions. Similar ratios can be taken for class II amplitudes. Based on the results of our 
model, we expect 

BR(B° -> KJ/ib) r 

^= mw^wim KMi[M2 ' °- 321 ' (32) 

The corresponding experimental value is R§ = 0.58 ± 0.11. 



3.2 B c decays 

The transition operators driving B c decays are the same that generate B and D decays. However 
their expectation values are evaluated for the B c wave function, rather than the B and D wave 
functions reflecting that b — > c and c — > s transitions proceed in a different environment. This 
is illustrated by the results of Table 4, where we compare various B c form factors at q 2 = 
with the corresponding B° and D form factors. 

We apply the strategy outlined in the previous sections to calculate different partial rates 
and the lifetime of the B c . The critical point with regards to this issue is charm decay. Here the 
energy release is not as comfortably large as it is in the case of bottom decay. As a result, our 
estimations of the inclusive charm decay should be more sensitive to a hadronization model. 
However, the inclusive charm decay contributes only 8% to the total c — > s rate of B c . For this 
reason we do not include any hadronization corrections in our calculations. 

The results for the partial B c decay modes corresponding to the various underlying quark 
subprocesses are collected in Table 5 for M% = Mp + m n . As in Table 2 we present the 
numbers of the exclusive and inclusive channels included in the calculations. The number of 
the exclusive b — > c decays in the B c case is less than that for B° because of the interference 
effects in the class III amplitudes (see Table 7)Q For comparison we also show in Table 5 
the results obtained in [jllj, using the nonrelativistic QCD. Viewing this comparison with due 



caution, regarding the model dependence and other uncertainties in the estimation of the decay 
modes as well as the quark mass uncertainty for the inclusive prediction, it is reassuring that the 
order of magnitude comes out to be consistent. Our bound state corrections are numerically a 
bit larger than very small effects reported in [fTOfl . For the sum of b — > c spectator contributions 
we obtain r^ c ) = 0.83 ■ T(B°) = 0.531 ps' 1 , while for the total c-decay contribution one 
finds r( c "^ = 0.38 • T(D°) = 0.915 ps" 1 . For the choice M^ 0) = M v + m n we obtain = 
0.79 • T(B°) = 0.504PS" 1 and T^ 3 ^ = 0.43 ■ T(D°) = 1.049 ps' 1 Q. Various branching fractions 
can be also inferred from Table 5. For instance the SL branching ratio BR(B C — > euX) is found 
to be 9.3%, i.e. ~ 15% less than BR(B° — > evX). Details of our predictions for the partial 
widths are presented in Table 6. In particular, we obtain 

T(B C -> ev e X- cc ) = 4.1 ■ 10 13 |K6| 2 sec-\ (33) 

Y(B C -»■ ev e X ls ) = 8.5 ■ l® lQ \V cs \ 2 sec~ l . (34) 



6 Note that in the heavy-quark limit the fraction ratios Ri, i = 1, ...4 equal unity. 

7 The sum of the other interference effects is effectively accounted for by the Pauli interference contribution 
Tpi, see below 

8 The b — > c rates have been calculated usig the values of the effective CKM parameter \V c b\ listed in Table 
2 for cases a and b, respectively 
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The corresponding results of the original ISGW2 model are [19 



T(B C -> ev e X- cc ) = 3.36 • 10 13 | V^sec' 1 , (35) 
Y(B C -> ei/ e X 5s ) = 5.0 • 10 10 |y cs | 2 sec -1 . (36) 
The ratio of B c decay via b — > c to c — ► s decay is 

r(B -> evA) = 0Q21 1M M L3 (37) 

T(5 C -> ez/ e X gc ) i h c i 2 

The SL width r(5 c -> ez/ e X 5c ) = 0.061 ps^ 1 is only » 20% smaller than that of 5°, 1(5° -> 
e^ e X g rf) = 0.074 ps -1 , even though the spectator c quark is no longer light. The contributions 
from the pseudoscalar and vector final states are 21% and 44%, that is not too much different 
from the predictions of the ISGW2 model. For the SL c — > s decays the recoil effects are very 
small due to the large daughter mass. The SL width T(B C — > eu e Xi s ) = 0.081 ps' 1 is dominated 
by the decays to the B Q S (27.2%) and B°* (63.2%), because the avaiable energy is small. 

In our phenomenological description we also include non-spectator contributions from weak 
annihilation (WA) and Pauli interference (PI) ||36|| . The the contribution of the annihilation 
channel is 

r a =£— \v a \*m* f^Vf^Vfi-A 2 V-Q, (38) 

where f Bc = 0.42 Gev, c T — 1 for the r + v T channel and c = (2c + (2/i re d) + c_(2/i re rf)) 2 /3 for the 
cs channel, with \x re <i = m^m c j [m^ + m c ) being the reduced mass of the bc system. The result 



is known to two-loop accuracy |37]: c + (2fi red ) = 0.8, c_(2/i rerf ) = 1.5. The leptonic constant 
fs c has been estimated through the LF wave function ip(x,p 2 ± ) using Eq. (13) of ref. 



Note that because of the partial cancellation of weak annihilation rate and the effect of the 
Pauli interference diagrams, no significant uncertainty on the lifetime arises from the limited 
knowledge of the decay constant f Bc . We find T a = 0.189 ps" 1 and T PI = —0.100 ps~ l . The 
similar value of Ypi has been obtained in [|l(|. Putting everything together we obtain for the 
B c width T(B C ) = T b -*' C (B C ) + r°-»(B e ) + T a + T PI : 

r(B c ) = (0.65 ps)-\ case a; F(B C ) = (0.61 ps)-\ case b. (39) 

One observes a dominance of the charm decay modes over 6-quark decays. We consider the 
dispersion in predicted values of t Bc as a rough measure of our theoretical uncertainty in 
calculation of the inclusive decay rates, therefore our final result is 

r Bc = 0.63 ± 0.02 ps~ l . (40) 

Note that the exclusive b- and c-decay rates are 0.32 ps" 1 and 0.80 ps -1 , i.e. 52% and 87% 
of the total b — > c and c — > s decay rates, respectively, In Tables 7,8 we compare our results 
for the exclusive two-body NL B c decay modes with the results obtained using the BSW and 
ISGW models |§ and with the results of ref. 0. Adding the rates for twenty of these two-body 
decay modes one finds 0.69 ps^ 1 (BSW model) and 0.90 ps^ 1 (ISGW model) and 1.15 ps^ 1 
in [0, to be compared with our result 0.70 ps" 1 which is the sum of the rates of 0.67 for 
c —>■ s transitions and 0.03 ps^ 1 for b — > c transitions from Table 7,8. 

Finally, we note that the experimental extraction of B c signal from the hadronic background 
requires the reliable estimation of the branching fraction BR(B C — > J/tp + X), because J/i[) 
can be easily identified by its leptonic decay mode, while the experimental registration of the 
final states containing the r/ c or B^ is impeded by the large hadron background. We obtain 
BR(B C -> J/ip + fiVp) = 1.7%, BR(B C -> J/^ + vr) = 0.1%, and BR(B C -> J/ip + X) = 13.2%. 



12 



4 Summary 



In conclusion, we have used a relativistic constituent quark model based on the LF formalism 
to perform a detailed investigation of the B c meson partial widths and the lifetime of the 
B c meson. The hadronic form factors and the distribution function were calculated using 
meson wave functions derived from an effective qq interaction intended originally to describe 
the meson mass spectra. In this way the link between B c physics and the "spectroscopic" 
constituent quark models was explicitly established. For numerical estimates we employed the 
LF quark functions that are related to the equal-time wave functions of the ISGW2 model. 
In several important aspects our analysis goes beyond the quark model estimations derived 
previously in the literature. In addition to the frequently discussed SL and the two-meson 
exclusive decay modes, we have included 84 exclusive decay modes and 44 SL and NL inclusive 
decay modes, that increases the total B c width by ~ 30%. Account of the inclusive decays leads 
to a considerably smaller £? c -lifetime, tb c = 0.63 ± 0.02 ps, which is in fairly good agreement 
with the estimate obtained non-relativistic QCD, r Bc = (0.55 ± 0.15) ps [|l(J and also agrees 
with the most recent CDF measurement H] within one standard deviation. 

To sum up, the LF constituent quark model makes clear predictions on the global pattern: 
(i) a short B c lifetime well below 1 ps and (ii) a predominance of charm over beauty decays 
among the NL modes. 

Appendix 

The vector and axial hadronic form factors relevant for weak SL decays of a pseudoscalar meson 
Hq, containing a heavy quark Q, to a member of the lowest-lying multiplet of pseudoscalar or 
vector mesons are defined in a usual way (see e.g. |30]]). We denote by Ph q , Ph q , and Mh q , 



Mh qI the 4-momenta and masses of the parent and daughter meson, respectively. Then the 
amplitude < Ph q ,\Vq\Phq > can be expressed in terms of two form factors 



Ml -Ml Ml -Ml 

< Ph q ,K\Ph q >= (P a - Q q2 Q q a )Fi(q 2 ) + Q q2 Q q*F (q 2 ), (A.1) 



where P = Ph q + Ph q , , q = Ph q — Ph qI ■ There is one form factor for the amplitude 

< Ph q ,, £ \Voi\Phq > 



< PHoii £ \^o\PHo > 



2i 



Mh q (1 + () 



'—-^se^PlptViq 2 ), (A.2) 



where ( = M Q f\ and three independent form factors for the amplitude < Pjj ,, e\A a \Pn Q > 



< Ph qI , £\A a \PH Q > 



M HQ (l + C)e^A 1 (q 2 
e*q) 



(e*q) 



2M, 



Mh q + M Hq 
e*q 



-P a A 2 (q 2 



q 



q Q A 3 (q 2 )+2M HQ -^A (q 



(A.3) 



9 In what follows we suppress the indeces 1 and 2 (corresponding to the meson containing a quark-spectator 
and to the one generated from a W current, respectively) of £ everywhere except in Eq. (A. 17). 
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where Aq(0) = ^3(0) and the form factor A^(q 2 ) is given by the linear combination 

Mi 2 ) = ^ ((1 + OMq 2 ) - (i - QA 2 { q 2 )) . (A.4) 

We view the form factors as functions of the dimensionless variable y = vh q • v H Q n where 
P Hq = M Hq ■ v Hq , P Hq , = M HqI ■ v HqI , and q 2 = M 2 Hq + M 2 Hqi - 2M Uq M HqiV Q. The 
differential rates for a decay of a pseudoscalar particle into another pseudoscalar particle and 
lepton or quark pair is given by 

L6 

-T c 2 \V QIQ \ 2 \V qiq2 \ 2 C%{q 2 ,mlml)J ■ 

V^IP - (V)(y 2 - 1) + ^ 2 (C 2 - (y)] (A.5) 

where Vq/q and V qiq2 are the corresponding CKM matrix elements and the factor J takes into 
account the perturbative QCD corrections. For the SL decays c 2 = 1, V qiq2 = 1, and J ~ 0.9; 
for the external (internal) NL decays c 2 = 3c 2 (3c 2 ,) and J = 1. The differential transition rate 
to a vector particle is 




' /J = 4r c 2 |^Q| 2 |^ l92 | 2 C 2 $(g 2 ,m 2 ,m 2 )J 



where 



Z/ 2 - 1[(1 - y)#%) + 6A a CV - l)A 2 (y)\, (A.6) 



# 2 (l/) = T^(H\{y) + + tf %)), (A.7) 



and the helicity amplitudes are given by 



c 



H ± {y) = M Hq \(l + QA^y) =F 2^ _ ly^VWJ , (A.8) 
M 2 ( C \ 

My) = Uy - C)(i + 0AM - 2( y 2 - i)J^a 2 (v)\ (a.9) 

The maximal value of y is (M# + Mfj )/(2M H a M H ), corresponding to g 2 = 0. 

Recall that in the case of the heavy-to-heavy transitions //q — > ifg/ the heavy quark 
symmetry implies simple relations between various form factors: 

Fi{y) = v{ y ) = My) = ^wMv), ( A -!0) 
F (y) = AM = ^ ( y+ 2 h iw(y), (A.11) 

where Ciw(y) is the Isgur-Wise function. For the realistic case of finite heavy-quark masses, 
these relations are modified by corrections that break heavy quark symmetry. These corrections 
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Although we are using the variable v\ ■ vi we are not treating the daughter quark as heavy. 
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have been analyzed for B° and D decays using the LF technique [P0[] , pifl . We use the same 
technique to calculate various form factors for the B c decays. The values of these form factors 
at q 2 = have been already given in Table 4. 

The NL two-body partial widths are given by expressions 

TftP, = y ■ T \V QV \ 2 \V qiq2 \ 2 F 2 PS ■ (1 - C 2 ) 2 «F 2 (C 2 ), (A.12) 

= y •r |^ Q | 2 |^ 1(?2 | 2 ^.« 3 F 1 2 (c 2 ), (A.i3) 

IftVi = j ■ UV Q/Q \ 2 \V qig2 \ 2 F 2 PS ■ a 3 A 2 (C), (AAA) 

= j ■ r |^ Q | 2 |K iq2 | 2 F 2 • aH 2 (( 2 ), (A.15) 

where 

a = «(Ci, CD = VU + Cf-Cf) 2 -^ 2 , (A.17) 

and c — ai for the external two-body decays, c = a 2 for the internal two-body decays. 

Finally, the partial widths of inclusive NL decays in which the final state contains a charged 
or neutral meson directly generated by a color-singlet current are given by 

d^px r c 2 F| 2 2 

— - 1 — ^ — I V^Q'Ql \V qi q 2 \ ■ 

[(1 + s- ( 2 ) 2 W 2 - 4C 2 W, + 4C 2 (1 + s - C 2 )W>(C 2 , s) (A.18) 
r 2 w 2 

^ = r ^|^ Q | 2 |^ ig2 | 2 -[a 2 (C 2 , S )^ 2 + 12C 2 ^]a(C 2 , S ) (A.19) 
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Table 1. The oscillator parameters (5 for various mesons |19| and the meson constants fps, fv 
[|30| used in the calculations. 



Meson 


IT 


K 


Vc 


D 


D s 


B° 


B s 


B c 


P 


K* 




D* 






0.41 


0.44 


0.88 


0.45 


0.56 


0.43 


0.54 


0.92 


0.30 


0.33 


0.62 


0.38 


0.44 


fps 


0.13 


0.16 


0.38 


0.21 


0.23 






0.42 












fv 


















0.20 


0.21 


0.41 


0.20 


0.27 



Table 2. The inclusive partial widths of the B° meson for the different choices of the contin- 
uum threshold M% (in units | Vf>c| 2 /|0.039| 2 ps^ 1 ). We use the values of the CKM parameters 
| Vac | — 0.974, \Vdc\ = \V US \ = 0.221. The CKM matrix element \V c b\ is calculated using the ex- 
perimental value of T tot (B°) = 0.641 ps~ x . Case a: Aff = My + m n , case b: Mjp = M P + m n . 



Decay mode 


N exc 


A^inc 


a 


b 


b —> c + e + v e 


2 


1 


0.074 


0.082 


b — ► c + ii + 


2 


1 


0.074 


0.081 


6 — > c + r + v T 


2 


1 


0.015 


0.017 


b — > c + w + J 


12 


6 


0.313 


0.358 


6 — > c + C + 5 


12 


6 


0.123 


0.135 


6 —> c + + s 


12 


6 


0.016 


0.018 


b —> c + c + J 


12 


6 


0.006 


0.006 


£° -> iVA c , A C S C 




2 


0.023 


0.023 


6 — > w 


54 


29 


0.007 


0.008 


^tot 


108 


58 


0.652 


0.728 


v bc 






0.0386 


0.0366 


Br S L 






11.35% 


11.26% 


n c 






1.198 


1.194 
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Table 3. The two-meson branching fractions of B° (in %). 



Decay mode 


This work 


ref.g3 


ref. m 


Exp. data 34 


Class I 










B° -> D + 7T~ 


0.34a? 


0.257a? 


318a? 


0.29 ±0.05 




0.78a? 


0.643a? 


0.778a? 


0.81 ±0.17 


B° L>*+7r- 


0.32a? 


0.247a? 


0.296a? 


0.25 ±0.05 


£° -> 


0.83a? 


0.727a? 


0.870a? 


0.70 ±0.16 




1.05a? 


0.879a? 


1.004a? 


0.82 ±0.36 


£° -> D+D*- 


0.92a? 


0.817a? 


0.830a? 


0.95 ±0.45 


B° D* + D; 


0.66a? 


0.597a? 


0.603a? 


0.85 ±0.33 


B° D* + D*~ 


2.32a? 


2.097a? 


2.414a? 


1.85 ±0.72 


B° -> D + K~ 


0.028a? 


0.020a? 


0.037a? 






0.045a? 


0.035a? 


0.041a? 




B u — > D* + K~ 


0.015a? 


0.019a? 


0.022a? 




B" —> D* + K*~ 


0.033a? 


0.042a? 


0.049a? 




Class II 










B° 7T°D° 


0.28a^ 


0.164a^ 


0.084a^ 


< 0.033 


B° 7T D 


0.21a^ 


0.230a| 


0.116a| 


< 0.055 


B° p°D° 


0.12a^ 


O.llla^ 


0.078a^ 


< 0.055 


B u — > p u Z}* u 


0.21a2 


0.240a2 


0.199a2 


< 0.117 


B° -> K° + rj c 


0.81a^ 








B° if J/V 


0.64a^ 


2.262a^ 


0.800a^ 


0.075 ±0.021 




0.41a^ 










1.46a| 


3.645a^ 


2.518a^ 


0.151 ±0.091 
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Table 4. Form factors for the external _B+-decays at q 2 = 0. In parentheses are given the 
values of corresponding form factors for the decays of B° or D°. 





B c -> ?7 C , J/^ 


B c — > -B s , 5* 


fi c -> D°, D°* 


B c ^ B°, B°* 


Fi(0) 


0.622(0.683) 


0.564(0.780) 


0.089(0.293) 


0.362(0.681) 


A o (0) 


0.477(0.678) 


0.447(0.730) 


0.050(0.214) 


0.282(0.600) 


A 1 (0) 


0.447(0.623) 


0.453(0.633) 


0.045(0.170) 


0.290(0.502) 


A 2 (0) 


0.390(0.556) 


0.522(0.464) 


0.040(0.155) 


0.370(0.366) 


V(0) 


0.621(0.683) 


2.17(0.777) 


0.073(0.216) 


1.50(0.663) 



Table 5. Inclusive B c partial rates (in units of ps 1 ). 



Decay mode 


N exc 


A^inc 


This work 


ref. [TO 


b — > c + e + v e 


2 


1 


0.061 


0.075 


b -> c + fi + 


2 


1 


0.061 


0.075 


b — > c + r + z/ r 


2 


1 


0.013 


0.018 


6 — > c + w + d 


12 


6 


0.259 


0.310 


b — > c + c + s 


8 


6 


0.102 


0.137 


6 — > c + ii + s 


12 


6 


0.013 


— 


b — > c + c + J 


8 


6 


0.006 


— 






2 


0.011 




o — > it 


46 


29 


0.005 




^tot 


92 


58 


0.531 


0.615 


c — ► s + e + z/ e 


2 


1 


0.081 


0.162 


c -> S + /J + ^ 


2 


1 


0.077 


0.162 


c — > s + w + J 


12 


6 


0.698 


0.905 


C — )> S + IX + s 


12 


6 


0.023 




c —> d 


28 


14 


0.036 




yc 
^tot 


56 


28 


0.915 


1.229 


be —> T + V T 






0.052 


0.056 


be —> c + s 






0.137 


0.138 


PI 






« -0.100 


-0.124 


^tot 


148 


86 


1.535 


1.914 
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Table 6. SL and NL external partial rates for b — > c and c — > s transitions (in units 10~ 3 ps _1 ). 
The values of the CKM parameters used in the calculations are the same as in Table 2. The 
values of the QCD two-meson amplitudes are a\ = 1.0, a 2 = —0.3 and a\ = Ci(m c ), a 2 = c 2 (m c ) 
for the b — > c and c — > s transitions, respectively. 





6 — > c/z// 


c — > s/z// 






X C c 


ft 


ft* 




ez/ e 


13.05 


26.6 


20.6 


22.0 


51.2 


7.35 




13.02 


26.5 


20.6 


21.1 


48.2 


6.34 


rz/ r 


4.37 


7.53 


1.81 











6 — > cud 


b — > cus 




P + 


X u d 


K+ 








2.25 


5.15 


46.4 


0.178 


0.289 


2.30 




1.27 


3.57 


95.9 


0.096 


0.213 


4.61 


x c - c 


5.71 


13.3 


65.9 


0.329 


0.751 


2.99 






b — > ccs 


b —> ccd 


ft 


ft* 




D+ 


L>*+ 


X c d 


Vc 


8.03 


6.68 


10.8 


0.344 


0.198 


0.690 


J/lp 


3.28 


13.2 


17.5 


0.146 


0.354 


1.16 




6.79 


14.0 


2.68 


0.426 


0.440 


0.245 






c — > swd 


C — > SSI/- 




P + 


X u d 




#*+ 




ft 


96.4 


65.4 


0.180 


8.06 


0.432 




ft* 


55.0 


340 




3.22 


0.213 






54.5 






1.56 
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Table 7. Two-meson b — > c rates of (in 10 3 ps 1 ). 



nppQv monp 

J—/ V_y VvCAj y 111 v./ VJ.V 


This work 


rpf ifil 


rpf l7l 
II Ml 


PlflQQ T 

Vyldoo 1 










o oq^2 

z.zocii 




Q 1 A ^ 2 

0. 14a^ 


D+ v r, n + 

&c ~^ VcP 




4 ^9n 2 




D+ , T/a/,vr+ 


1 9^/7 2 


1 Q9/7 2 




D+ > T/ohn+ 

&c ~^ J /VP 






Q (14/7 2 


£> c — ► ?7 C A 


n oq^2 

U.ZoO^ 


n 1 /i /i 2 

U. 14(2-^ 


n 9/i e;^ 2 

U.Z40(2^ 


JD C — > A/ c 7\ 


3(S<7 2 

U.UUUii 


n 99o 2 


D 4S C: )(7 2 




12(7? 

17 . _L Z- 1.4/ 1 


14(7? 

V 7 . -L Til. 4 1 


231(7? 


p>+ _^ J/ihK* + 


27a? 

V 7 . Z- 1 Liy 1 


28a? 


492a? 

V 7 • x i_/ Zj La/ 1 


PlflQQ TT 








ij c — >■ u u 


/I 1 Q^ 2 
4. L6CL2 




1 m ^ 2 
14JICI2 




3 1Qo 2 




1 D(So 2 
1 . UUU2 




1 31fln 




992«? 




2 01(7? 




1 64a? 


PIqqq ttt 

Lyldbb 111 










/q 1 n~ 1 /i 1 n~ \2 

^o.iyai + 4.iya2j 




(^i.oyai 4- z. 44G,2j 




f9 Q1 j_ 9 4^/7^ 2 




yL.Zoui ~t~ Z.o4tt2j 


D+ v T /oh 74+ 






^l.ZDtii T" Z.4Uti2j 


d+ v 7/,/, n*+ 


^4.Uytti ~r 4. / 0U2 ) 






r> c — > r/ c i^ 






rn 94/^1 -i- n ^dn n \ 2 




0.45ai 4- 0.81a 2 ) 2 




(0.22ai + 0.53a 2 ) 2 


Bt - 


(0.37ai + 0.54a 2 ) 2 




(0.22ai + 0.54a 2 ) 2 




(0.59ai + 0.98a 2 ) 2 
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Table 8. Two-meson c — > s rates of (in 10 3 ps 1 ). 



Opppv TVTnrlp 

L/V A ( l y 1 V J_ v 7 V J. v 


This work 


rpf ifil 


rpf ifil 


lei. m 


PlflQQ T 

Vyldoo 1 










i^T — > # S 7T^ 


CO 0£? ^2 

oz.oba 1 


A 1 o^2 

4/.za 1 


cc o ^2 

ob.oa 1 


oo 7^2 

oo. /a| 


#c B sP 


oo.ooa 1 


1 n n ^»2 

iy.Ua 1 


on 7^,2 

oU. /a x 


<^o 1 ^»2 

oo.la 1 




on ~\ a ^2 
o(J.14a 1 




rq 7^2 

oz. / a x 


70 /| .,2 

/o.4a 1 




1 /i a J2 

lob.4oa 1 


1 /o.ba 1 


oo 1 1^2 

zol.la 1 


ooo „2 

ZZ6a 1 


D c — > UqTI 


o qn~2 

z.oua^ 


i /| 7^2 
1.4 / a^ 


z.o < a^ 


^ m ^2 
o.uia^ 




Zj. JUUiI 


1 4S<7 2 




Q D7o 2 




1 19(7? 


2 40tf? 


1 Q4<7? 


4 41(7? 




10 3a? 

IV/. ' /V A | 


13 4a? 


13 5a? 


18 07a? 


PIoqq TT 










ii c — > U K 


of? e /( ~2 

oo.o4a 2 


/i o o ^,2 

4z.oa 2 


no q^,2 

yo.oa 2 


i /i n J2 

14:1 a 2 


— > UK 


on 

zU.9oa 2 


i r O ^2 

15. za 2 


ob.ba 2 


i nj „2 

lU4a 2 


D+ > D*+ Z^O 

UJ — > U K 


in rr„2 

lo.ooa 2 


A 7 1 „2 

4/.la 2 


AO C\ ^2 

4o.Ua 2 


111 „2 

llla 2 


li r T — > U K 


inr i n„2 

lzo.lUa 2 


ooo c ^2 

zzo\oa 2 


o a o n ^2 

z4o.ya 2 


01 A sJ2 

z!4a 2 


ft + — > H + 7T 

C 


1 . <JUCt/ 2 


n 7^ 2 


1 44^ 2 
1 . z ± z ±(x 2 


9 ^1 n 2 

Z. U lCt' 2 




1.95a 2 , 


0.71a^ 


1.63a 2 , 


4.53a^ 




0.80a 2 , 


1.20a 2 , 


0.97a^ 


2.20a^ 


5+ -> 


6.93a 2 , 


6.70a 2 , 


6.73a^ 


9.05a^ 
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Figure 1: Distribution functions -Fg(x) for B° and B c mesons. 
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